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The influence of earth curvature in the theory of reflection from the ionosphere is con- 
sidered. By choosing a rather idealized model, the significance of usual earth-flattening 
procedures can be displayed quite readily. To simplify the analysis, it is assumed that the 
earth's magnetic field is vertical everywhere. It is shown that the curved ionosphere may be 
represented by homogeneous planar slabs, provided the local value of layer curvature is used. 
The results in the present paper are compared with some corresponding expressions obtained 
by Johler and Berry [1962]. 



1. Introduction 

Boundary value problems involving magneto 
plasmas, such as the ionosphere, are a great deal more 
complicated than the corresponding problems for 
isotropic media. For this reason a number of gross 
simplifications are often made in order to achieve 
tractability. For example, in most of the existing 
literature on radio propagation it is assumed that the 
ionosphere may be regarded as a planar stratified 
medium. In fact, Budden's exhaustive treatise 
[Budden, 1961] is restricted entirely to this case. 
Also, in prior work on VLF mode theory it has been 
somewhat tacitly assumed that the ionosphere may 
be regarded, in a local sense, as a planar stratified 

.medium. Justifications for this step are based on 
arguments valid only for isotropic media. Earth 
curvature is accounted for only in the nonionized 
waveguide region beneath the base of the ionosphere 
[Wait, 1961a, 1962; Katzin and Koo, 1962; Budden, 
1962]. Since the anisotropic properties of the iono- 
sphere play an important role in VLF mode theory, it 
seems that this aspect of the subject warrants some 
scrutiny, and this is the main purpose of the present 
paper. 

The situation considered here is rather idealized. 
The ionosphere is regarded as a spherically stratified 
plasma with a superimposed radial magnetic field of 
constant strength. Such a model has been considered 
by van der Wijck [1946], Bremmer [1949], and more 

^recently by Krasnushkin [1961]. If this restriction, 
on the magnetic field, were not made the problem 
would become intractable unless some other simplify- 
ing assumptions were made such as a purely trans- 
verse field [Wait, 1961b]. 

2. Formulation 

Choosing a spherical coordinate system (0, 0, r), 
the base of the ionosphere can be regarded as the 



spherical surface r=c. For an ionospheric-type 
plasma, the dielectric constant (e) has the form of a 
tensor. The relation between the displacement 



vector D and the electric field E is then 
D=(e)E, 
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which may be written more explicitly in matrix form 
as 



(2) 



The elements of the dielectric tensor may be a 
function of r but constant with respect to 6 and 0. 
They depend on the strength of the earth's magnetic 
field and on the electron density, ion density, and 
frequency of collisions between them. The case 
usually considered is that in which the electromag- 
netic forces influence only the electrons and the 
motion of the ions is neglected. In this case 
[Bremmer, 1949] 
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where w is the (electron) plasma frequency, w h is 
the (electron) gyrofrequency, and v is the effective 
collision frequency (for electrons). The generaliza- 
tions of these formulae to include the influence of 
heavy ions and the effect of an energy-dependent 
collision frequencv have been discussed elsewhere 
[Sen and Wyler/l960; Johler, 1962; Wait, 19621. 
For purposes of the present analysis, it is sufficient 
to point out that, even under these general conditions, 
the dielectric tensor has the form given by (2), 
provided the magnetic field is vertical. Therefore, 
the subsequent theory will not be restricted just to 
an electron plasma with constant collision frequency. 
Maxwell's equations, 



curl E=—ifjiajH, 
curl H=-i(e)o)E, 



(4) 
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where (e) is a tensor may be written explicitly in the 
form 
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In the above it has been assumed that d/dc£=0. 

It is now a simple matter to express the r and 6 
components of the field in terms of the <f> component. 
To facilitate this, two auxiliary functions M and TV 
are introduced which are defined by 
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Furthermore, it is found that M and A T satisfy the 
two coupled equations 
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(sin0|M)= 



= 0. (19) 



where k 2 = i nw i . In the case of isotropic media, 
M and N are proportional to Debye potentials. 
Equations (18) and (19) are applicable in regions 
external to the source. 

3. Solution 

A solution of (18) and (19) is now sought in the 
form 

N=F(r)P v (-cos6), (20) 



< 



M=G{r)P v {- cos- 
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where P v is the Legendre function which has the 
required asymptotic behavior and has no singularities 
except at 0=0. It is a solution of the equation 
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The radial functions F and G must then satisfv 
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In the case of isotropic media (i.e., e 12 = 0), the equa- 
tions become uncoupled. As a further check, the 
medium may be taken to be homogeneous and iso- 
tropic (i.e., eii=€22=€ 3 3=€ where e is a constant). 
Then F and G satisfy the same differential equation 
given by 
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Solutions are linear combinations of the functions 
Zh£\Z) and ZK 2) (Z), where Z=k(e/eo) 1/2 r and where 
/>i n and h [ v 2) are spherical Hankel functions. Most 
proponents [Budden, 1962; Katzin and Koo, 1962] 
of earth-flattening procedures start with an equation 
equivalent to (25). Then they develop the idea that 
the term 



k 2 



may be replaced by 



e 



■m 



I where e eff (r) is the effective dielectric constant for 
: the equivalent planar region. Unfortunately, in the 
anisotropic case the solutions of the coupled equations 
cannot be expressed in terms of spherical Hankel 
functions even if the medium is homogeneous. 
Thus, the earth flattening procedures are open to 
question when dealing with magnetoplasma media. 
To solve the coupled wave equations for the genera] 
case, resort must be made to numerical techniques. 
To simplify the procedure it is suggested that the 
ionosphere above the surface r=c=r be broken into 
a number of concentric homogeneous regions which 
are bounded by the surfaces r , r u r 2 . . . r m -i, r m . . . 
r M . For example, within the rath region, the ele- 
p ments of the dielectric tensor may be replaced bv 
constants. Thus, for T m - X <ir<^r mt 

eii = e 2 2=e m , e 33 = € w and ei2=ig m , 

which corresponds to a dielectric tensor having the 
form 
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for the rath region. The coupled equations now 
have the form 
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a constant within the 4 layer. Thus, it is replaced by 
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for the rath layer. Solutions of the coupled differ- 
ential equations may now be found in the form 



and 
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G m =q m e-v, 
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Substituting these into (27) and (28) leads to the 
two algebraic equations 
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In order for a solution to exist, the determinant of 
the coefficients p m and q m must vanish. This 
condition can be written compactly in the form 

utn+uilsl+sll+slst-sl^^O. (34) 

4=^r(lr e '"-^> (35) 



where 



«?= 



-k 2 



•XL 



(36) 



There are four solutions of this equation; two 
representing outgoing waves and two representing 
incoming waves. These are denoted ±u m and ±u' m . 
The equation is equivalent to the Booker quartic for 
plane waves in anisotropic media [Budden, 1961]. 

The general solution in the mth layer may then be 
written 



and 



F m =p m e'^ r +p m e + ^ r +j^e'^ r +pUe + ^ r (37) 
^=ie-^ r +J m 6 + v+J;e-^ r +£e +u ^, (38) 



when the p's are arbitrary constants. The q's are 
related to corresponding p's through (32) or (33). 
Therefore, for convenience in what follows, 

-^ -> <- -> 

q m = a( U m ) p m , q m =<*( — U m ) p m , 

q'm=a(u m )p m , q' m =a( — Un)qm, 

where the function a(u) is known. The imaginary 
part of the u^s is chosen to be positive and, thus, 
the arrows -> and <- represent outgoing and ingoing 
waves, respectively. In the outermost region, where 
ra=M,only outgoing waves are permissible and, thus, 

<r~ <- <- «- 

p M = p' M =q M =qM~0. Since the tangential field com- 
ponents are to be continuous at the interfaces, it 
follows, from (12), (13), (16), and (17) that four 
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continuity conditions are required at the surface 
r=r m . These are 

F =F _i_i 

x m ■*- m-f-lj 



to m to m +l) 

(lG m dG m +i 



dr 



dr 



(39) 
(40) 

(41) 



Therefore, four independent relations between the 
coefficients are available at each interface from 
m=l to m=M. 

To obtain the complete solution of the problem it 
is necessary to consider the fields in the free space 
or waveguide region a<^r<^r . Here the fields are 
obtained from 
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which is a special case of equations (12) to (17) for 
isotropic media. In this case 



N =F (r)P,(- cos 0), 



and 



M = G (r)P,(- cos 
where G and F satisfy 
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Therefore, solutions in the region a^r^r^ are of 
the form 
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The boundary conditions at the level r=r a are 
F =F 1 
G = c?i, 
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For all practical purposes, the ground may be 
characterized by a surface admittance Y g and a 
surface impedance Z g . Thus 



and 
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at r=a. In terms of G and F these may be written 
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at r=a. In the case of a homogeneous ground of 
conductivitv a« and dielectric constant e„, 
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In most actual situations the square bracket factors 
in these expressions may be replaced by unity, and, 
furthermore, displacement currents in the ground 
are usually of minor significance. Thus 



Z g c±l/Y g e±(in<a/<r g y 
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The boundary conditions as exemplified bv (39), 
(40), (41), (42),* (60), and (61) yeild 4A/+2 algebraic 
equations connecting 4M+2 unknown coefficients. 
Setting the determinant of these coefficients equal 
to zero leads to an explicit equation to solve for the 
eigenvalues or modes p. To illustrate the form of 
this modal equation, the case of a homogeneous 
ionosphere is considered. Thus M=l and the 
boundary conditions yield 6 linear equations in the v 

coefficients /i, / 2 , <7i, 9-2, Pu and p[. The determinant 
of these coefficients is written out explicitly in table 
1. Setting this monster equal to zero yields a rather 
involved transcendental equation for the determina- 
tion of the modes v. The size of the determinant 
grows linearh- as the number of slabs is increased 
(i.e., M>1). In the case of any number of slabs, 
the determinant equation for the mode can be 
written in matrix form as follows: 
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Table 1. Determinant of the coefficients for the case of homogeneous ionosphere 
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Note: In this determinant M**> (Z)=Zh v ^ (Z), h v ^'{Z) =— [Z/i„(«) (Z)], where i=l, 5 
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where the first two factors arc matrix reflection 
coefficients. Explicitly, 
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where In' is the logarithmic derivative defined by 

dF(x)/dx 
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7?i and 74 can be interpreted as ground-reflection 
coefficients for vertically and horizontally polarized 
waves, respectively. The explicit form of \\R\\, 
j_R\\, etc., for the present problem could be written 
out in full. However, there is no need to do this 
as they have the same form as the reflection coef- 
ficient matrix for a plane stratified ionosphere. 
This aspect of the subject has been discussed recently 
in some detail [Wait, 1963]. To make direct use 
of numerical data already available for planar 
reflection coefficients, it is necessary to invoke the 
Debye approximations for the Hankel functions 
h ( p v (kc) and h { v 2> (kc) and their derivatives. The 
restriction for this step is 

where C is the (complex) cosine of the local angle 
of incidence at the ionosphere. In terms of v and 
kc it is not difficult to show that 
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In the theoretical model discussed in this paper, 
the earth's magnetic field was taken to be vertical 



everywhere. A rigorous treatment of the earth- 
ionosphere waveguide with a nonvertical magnetic 
field does not seem to be possible. However, an 
approximate approach to the subject, which utilizes 
the concept of surface impedances, also leads to a 
modal equation of the form of (65) above [Wait and 
Spies, 1960; Wait, 1961a, 1963]. 

4. Final Discussion 

The particular problem discussed in this paper has 
also been considered by Johler and Berry [1962]. 
However, their solution for the anisotropic ionosphere 
does not agree with the one given in this paper and 
with a previous formulation of the author [Wait, 
1961a] using a similar approach. It appeal's that 
the origin of the discrepancy is their assumption 
that the Hertz potentials IP and U m in the ionosphere 
are not coupled by the boundary conditions. 1 This 
coupling is evidenced by the forms of (18) and (19) 
in the present paper. Furthermore, their Hertz 
potentials IT and U m , which are related to the Debye 
potentials N and M, are not solutions of uncoupled 
wave equations in the ionosphere. 

The coupling of ordinary and extraordinary waves 
by the boundary conditions appears to be an essential 
point in the analysis. Physically this is important 
because the transmitted waves in the ionosphere are 
almost circularly polarized even though the waves 
in the space between the earth and the ionosphere 
are nearly linearly polarized. 
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